Whole-genome resequencing confirms reproductive isolation between sympatric demes of brown trout (Salmo trutta) detected with allozymes.

The sympatric existence of genetically distinguishable populations of the same species remains a puzzle in ecology. Coexisting salmonid fish populations are known from over 100 freshwater lakes. Most studies of sympatric populations have used limited numbers of genetic markers making it unclear if genetic divergence involves certain parts of the genome. We returned to the first reported case of salmonid sympatry, initially detected through contrasting homozygosity at a single allozyme locus (coding for lactate dehydrogenase A) in brown trout in the small Lakes Bunnersjöarna, Sweden. First, we verified the existence of the two coexisting demes using a 96-SNP fluidigm array. We then applied whole-genome resequencing of pooled DNA to explore genome-wide diversity within and between these demes; nucleotide diversity was higher in deme I than in deme II. Strong genetic divergence is observed with genome-wide FST  ≈ 0.2. Compared with data from populations of similar small lakes, this divergence is of similar magnitude as that between reproductively isolated populations. Individual whole-genome resequencing of two individuals per deme suggests higher inbreeding in deme II versus deme I, indicating different degree of isolation. We located two gene-copies for LDH-A and found divergence between demes in a regulatory section of one of these genes. However, we did not find a perfect fit between the sequence data and previous allozyme results, and this will require further research. Our data demonstrates genome-wide divergence governed mostly by genetic drift but also by diversifying selection in coexisting populations. This type of hidden biodiversity needs consideration in conservation management.

A Population Genomic Investigation of Immune Cell Diversity and Phagocytic Capacity in a Butterfly.

Insects rely on their innate immune system to successfully mediate complex interactions with their internal microbiota, as well as the microbes present in the environment. Given the variation in microbes across habitats, the challenges to respond to them are likely to result in local adaptations in the immune system. Here we focus upon phagocytosis, a mechanism by which pathogens and foreign particles are engulfed in order to be contained, killed, and processed. We investigated the phenotypic and genetic variation related to phagocytosis in two allopatric populations of the butterfly Pieris napi. Populations were found to differ in their hemocyte composition and overall phagocytic capability, driven by the increased phagocytic propensity of each cell type. Yet, genes annotated to phagocytosis showed no large genomic signal of divergence. However, a gene set enrichment analysis on significantly divergent genes identified loci involved in glutamine metabolism, which recently have been linked to immune cell differentiation in mammals. Together these results suggest that heritable variation in phagocytic capacity arises via a quantitative trait architecture with variation in genes affecting the activation and/or differentiation of phagocytic cells, suggesting them as potential candidate genes underlying these phenotypic differences.

Ecology and Genetic Structure of the Parasitoid Phobocampe confusa (Hymenoptera: Ichneumonidae) in Relation to Its Hosts, Aglais Species (Lepidoptera: Nymphalidae).

The biology of parasitoids in natural ecosystems remains very poorly studied, though they are key species for their functioning. Here we focused on Phobocampe confusa, a Nymphalini specialist, responsible for high mortality rates in charismatic butterfly species in Europe (genus Aglais). We studied its ecology and genetic structure in connection with those of its host butterflies in Sweden. To this aim, we gathered data from 428 P. confusa individuals reared from 6094 butterfly larvae (of A. urticae, A. io, and in two occasions of Araschnia levana) collected over two years (2017 and 2018) and across 19 sites distributed along a 500 km latitudinal gradient. We found that P. confusa is widely distributed along the latitudinal gradient. Its distribution seems constrained over time by the phenology of its hosts. The large variation in climatic conditions between sampling years explains the decrease in phenological overlap between P. confusa and its hosts in 2018 and the 33.5% decrease in the number of butterfly larvae infected. At least in this study, P. confusa seems to favour A. urticae as host. While it parasitized nests of A. urticae and A. io equally, the proportion of larvae parasitized is significantly higher for A. urticae. At the landscape scale, P. confusa is almost exclusively found in vegetated open land and near deciduous forests, whereas artificial habitats are negatively correlated with the likelihood of a nest to be parasitized. The genetic analyses on 89 adult P. confusa and 87 adult A. urticae using CO1 and AFLP markers reveal a low genetic diversity in P. confusa and a lack of genetic structure in both species, at the scale of our sampling. Further genetic studies using high-resolution genomics tools will be required to better understand the population genetic structure of P. confusa, its biotic interactions with its hosts, and ultimately the stability and the functioning of natural ecosystems.

Physiological Tradeoffs of Immune Response Differs by Infection Type in Pieris napi.

Understanding the tradeoffs that result from successful infection responses is central to understanding how life histories evolve. Gaining such insights, however, can be challenging, as they may be pathogen specific and confounded with experimental design. Here, we investigated whether infection from gram positive or negative bacteria results in different physiological tradeoffs, and whether these infections impact life history later in life (post-diapause development), in the butterfly Pieris napi. During the first 24 h after infection (3, 6, 12, and 24 h), after removing effects due to injection, larvae infected with Micrococcus luteus showed a strong suppression of all non-immunity related processes while several types of immune responses were upregulated. In contrast, this tradeoff between homeostasis and immune response was much less pronounced in Escherichia coli infections. These differences were also visible long after infection, via weight loss and slower development, as well as an increased mortality at higher infection levels during later stages of development. Individuals infected with M. luteus, compared to E. coli, had a higher mortality rate, and a lower pupal weight, developmental rate and adult weight. Further, males exhibited a more negative impact of infection than females. Thus, immune responses come at a cost even when the initial infection has been overcome, and these costs are likely to affect later life history parameters with fitness consequences.

Microevolutionary selection dynamics acting on immune genes of the green-veined white butterfly, Pieris napi.

Insects rely on their innate immune system to successfully mediate complex interactions with their microbiota, as well as the microbes present in the environment. Previous work has shown that components of the canonical immune gene repertoire evolve rapidly and have evolutionary characteristics originating from interactions with fast-evolving microorganisms. Although these interactions are likely to vary among populations, there is a poor understanding of the microevolutionary dynamics of immune genes, especially in non-Dipteran insects. Here, we use the full set of canonical insect immune genes to investigate microevolutionary dynamics acting on these genes between and among populations by comparing three allopatric populations of the green-veined white butterfly, Pieris napi (Linné; Lepidoptera, Pieridae). Immune genes showed increased genetic diversity compared to genes from the rest of the genome and various functional categories exhibited different types of signatures of selection, at different evolutionary scales, presenting a complex pattern of selection dynamics. Signatures of balancing selection were identified in 10 genes, and 17 genes appear to be under positive selection. Genes involved with the cellular arm of the immune response as well as the Toll pathway appear to be enriched among our outlier loci, regardless of functional category. This suggests that the targets of selection might focus upon an entire pathway, rather than functional subsets across pathways. Our microevolutionary results are similar to previously observed macroevolutionary patterns from diverse taxa, suggesting that either the immune system is robust to dramatic differences in life history and microbial communities, or that diverse microbes exert similar selection pressures.